Diamagnetism of doped two-leg ladders and probing the nature of their commensurate phases 



G. Roux,' El E. Origr 



and D. Poilblanc 



o 
o 

(N 
> 

o 
in 



a 

o 
o 



> 
o 

rn 
O 

o 



X 



' Laboratoire de Physique Theorique, IRSAMC, Universite Paul Sabatier, CNRS UMR5152, 
118 Route de Narbonne, 31400 Toulouse, France. 
^Laboratoire de Physique de I'Ecole Normale Superieure de Lyon, ENS-Lyon, CNRS UMR5672, 

46 Allee d'ltalie, 69007 Lyon, France. 
^Department of Physics and Astronomy, University of California, Irvine, California 92697, USA. 

(Dated: February 1, 2008) 

We study the magnetic orbital effect of a doped two-leg ladder in the presence of a magnetic field component 
perpendicular to the ladder plane. Combining both low-energy approach (bosonization) and numerical simu- 
lations (density-matrix renormalization group) on the strong coupling limit (t-J model), a rich phase diagram 
is established as a function of hole doping and magnetic flux. Above a critical flux, the spin gap is destroyed 
and a Luttinger liquid phase is stabilized. Above a second critical flux, a reentrance of the spin gap at high 
magnetic flux is found. Interestingly, the phase transitions are associated with a change of sign of the orbital 
susceptibility. Focusing on the small magnetic field regime, the spin-gapped superconducting phase is robust 
but immediately acquires algebraic transverse (i.e. along rungs) current correlations which are commensurate 
with the 4fcF density correlations. In addition, we have computed the zero-field orbital susceptibility for a large 
range of doping and interactions ratio J/t : we found strong anomalies at low J/t only in the vicinity of the 
commensurate fillings corresponding to 5 = 1/4 and 1/2. Furthermore, the behavior of the orbital susceptibil- 
ity reveals that the nature of these insulating phases is different: while for S = 1/4 a. Akp charge density wave 
is confirmed, the 5 = 1/2 phase is shown to be a bond order wave. 

PACS numbers: 74.20.Mn, 71.10.Pm, 75.40.Mg, 75.20.-s 



I. INTRODUCTION 

Ladder systems have proven to be remarkably interesting 
systems, both as simple models exhibiting behavior simi- 
lar to 2D systems, and as systems exhibiting competition 
between several types of ground states. Theoretical mod- 
els of doped ladders display a large superconducting (SC) 
pjj^g^ij2j3j4j5^ with d-wave pairing associated with the pres- 
ence of a spin gap, with a ground state which can be described 
variationally as a short-ranged resonating valence bond state. ^ 
Charge density wave (CDW) correlations compete with the 
pairing correlations j^i^i^ The phase diagram of the isotropic 
t-J model was sketched in Ref. 6 and displays, in addition to 
this competition, insulating phases for the particular commen- 
surate dopings 6=1/4 and 1/2. Another competition exists 
between the superconducting phase and an orbital antifeiTO- 
magnetic flux (OAF) phase^- '"''^''^'^^-'^ which has been ad- 
dressed in different ladder models by studying transverse cur- 
rent coiTelations which display a quasi-long range order in the 
OAF phase. 

Ladders are also among the simplest systems through which 
a magnetic flux can pass (see Fig. [T]l. When a magnetic field 
H is applied to an electronic system, it couples to both the 
spin of the electron, via Zeeman effect, and to the charge de- 
gree of freedom, via orbital effect. The total magnetic sus- 
ceptibility of a real material splits into various contributions'^ 
^ ^ X''''"+XTOnd+Xrare' whcrc x"'''" is Pauli susceptibiUty and 
Xcond Xcme ^'"^ Orbital susceptibilitics of the conduc- 
tion and the core electrons. Xcme must be evaluated from local 
atomic orbital and we will neglect it in what follows. x°md 
usually difficult to evaluate because one has to precisely de- 
scribe the evolution of the whole band structure with magnetic 
fieldfi^^ In the following, we will investigate Xcond within a 
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FIG. I : (Color online) The isotropic t-J ladder under magnetic field. 
If the magnetic field has a component perpendicular to the ladder 
plane, a flux (j) passes through each plaquette. Below is the gauge of 
the Peierls substitution with opposite phases ±<^/2 along the legs. 



single-orbital model. When the magnetic field is applied par- 
allel to the plane of the ladder, the orbital effect is suppressed 
and only the Zeeman effect remains. The latter case has been 
discussed in details in this system and it was shown that a 
doping-controlled magnetization plateau and a large Fulde- 
FeiTell-Larkin-Ovchinnikov phase are obtained. '^ '^•^'^ Since 
the ladder possesses a spin gap, when the magnetic field is not 
in the ladder plane, the orbital contribution may dominate the 
spin contribution in the total susceptibility at low temperature. 

Early numerical investigations of the t-J model with mag- 
netic orbital effect on ladders and 2D lattices revealed a strong 
effect of the magnetic field on the magnetic and pairing prop- 
erties,^' but the results were limited to small systems. A 
bosonization study of a related model of spinless fermionic 
ladders suggested the possibility of fractional excitations and 
of an OAF phase induced by the magnetic field. ^^'^^ Carr and 
Tsvelik"^ studied the orbital effect of the magnetic field on the 
interladder coupling using an effective model to describe a 
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single ladder. Lastly, it has been predicted^'' that bosonic lad- 
ders could have commensurate vortex phases at commensu- 
rate fluxes which would represent a one-dimensional analogue 
of the two-dimensional vortex phase. 

The purpose of the present paper is to consider the effect 
of a nonzero flux on the magnetic susceptibility on a single 
two leg ladder, and also to investigate the effect of stronger 
fluxes on the zero temperature phase diagram of the ladder. To 
this end, we combine the bosonization technique and density- 
matrix renormalization group^^'^^'^^ (DMRG, see AppendixlBl 
for details) to compute the phase diagram and physical prop- 
erties of a spin- 1/2 fermionic ladder with orbital effect. The 
results are presented in two parts. The first part is devoted to 
the analysis of the phase diagrams as a function of the flux 
(f) per plaquette in the weak- and strong-coupling limits. The 
second part is focused on the physics of the spin-gapped phase 
of doped ladder at small flux which is more relevant to real- 
istic magnetic fields. In particular, we discuss the stability 
of the insulating phases at (5 — 1/4 and 1/2 present in the 
phase diagram of the t-J model at zero flux. In the conclusion, 
we briefly give considerations on experiments which are con- 
nected to these results. For sake of clarity, we have relegated 
some of the technical details to appendices. 



Including the flux 

The magnetic flux couples to the kinetic part of the Hamil- 
tonian through Peierls substitution. In what follows, dif- 
ferent hopping amplitudes along the chains (<||) and between 
the chains (i^) are considered, and the Hamiltonian is: 



= -H J2 y'^^'4+l,l,aC^,l,a + h.C 
-h E [^~"^^^4+1.2,.^^.2,. + h.C 



(1) 



where Cij^a is the electron creation operator at site i on leg I 
with spin a. (f> denotes the dimensionless flux per plaquette 



^ * A{x)dx = -^H^a^ 



(2) 



with A{x) the vector potential which depends on the gauge 
choice, a the lattice spacing. The magnetic field breaks time- 
reversal and chain exchange symmetries which, as expected, 
will have notable consequences on the current properties of 
the system. Symmetry and periodicity considerations allow 
us to limit the study to flux < < tt. Exchanging chains 
amounts to reversing the direction of the magnetic field. More 
details on the gauge and flux quantization on a finite systems 
can be found in Appendix iBl The unit of </> is 27r0o with 
(f>o — h/e — 4.1357 x 10^^^ T m^. For experimental consid- 
erations, a flux ip — O.OItt akeady corresponds to a very high 
magnetic field of H 800 T for a typical value a = 4 A of 
the Cu-O-Cu bond in a cuprate. The gauge chosen in Eq. ([T]i 
is represented on Fig.[T] 



II. PHASE DIAGRAMS 



A. Weak-coupling limit 



We first introduce interactions between electrons using the 
Hubbard ladder in magnetic flux. The Hamiltonian comprises 
the kinetic term Ht incorporating the flux and the on-site re- 
pulsion U > 0: 



(3) 



According to the usual strategy,^ we will consider first the 
limit U — and study the non-interacting band structure. 
Then, we will turn on U ^ t^^,t± so that the band structure 
is not deformed and obtain the different phases in this weak- 
coupling limit. Let us begin with the discussion of the band 
structure. 

The magnetic flux has a strong effect on the shape of the 
bands. Indeed, it mixes the bonding (0) and antibonding (tt) 
bands which exist at zero flux. To emphasize the difference 
with the zero field case, we call the two bands obtained at 
finite flux the down (d) and up (u) bands. Results on the band 
structure are discussed in Refs. 22.23 and extended to a finite 
and fixed filling in Appendix lAl The band structure depends 
on the flux and on the ratio t±/t^^ (see Fig.|2]i. We define 
two characteristic fluxes (pc and 0o, both of them dependent on 
tj^/t||, and such that above (f>c, a double well appears in both 
bands and above (po, a band gap opens between the bands u 
and d. Four different possible shapes of the bands are obtained 
according to the location of the flux with respect to <j>c and (po- 
The two critical flux lines crosses at (t±/<|| — V^, 4> = 7r/2). 
By filling these bands, one can show (see Appendix lAl) that 
only situations with either two or four Fermi points can occur. 
The location of these Fermi points and their respective Fermi 
velocities vary continuously with the flux. In the rest of the 




FIG. 2: The four typical shapes of the bands depending on the flux 
and on the ratio tx/ty . Critical fields po (resp. pc) signal the appear- 
ance of a double-well (band gap). Note that the D phase always has 
only 2 Fermi points whatever the filling and that the C phase when 
p = n always has 4 Fermi points. 
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FIG. 3: Phase diagram in the weak-coupUng hmit for an isotropic 
ladder restricted to fillings < n < 1. Phases with 4 (resp. 2) Fermi 
points fall generically into the CISO (resp. ClSl) class. At half- 
filling, interactions will drive the system into a Mott insulating phase 
for < 00 while a band insulating phase occurs when (j> > (j)c. 
Other phases can be found if the ratio of the Fermi velocity is large 
(C2S1 and C2S2) and if the d-band Fermi wave- vector is 7r/2 (COSl 
andClS2). 



paper, we will work only at fixed electronic density (denoted 
by n) which will constrain the sum of the Fermi momenta as 
a result of the Luttinger theorem. ^"' 

Having obtained the non-interacting band structure, we add 
interactions small enough not to perturb the band structure, 
following the strategy of Refs. I4ll3lll32ll33ll3 4l Adopting the 
usual notation CpSq for a phase with p gapless charge modes 
and q gapless spin modes, a system with two Fermi points 
and repulsive interactions is expected to be generically in a 
ClSl phase, i.e. a Luttinger liquid state. With four Fermi 
points and repulsive interactions, the system is generically in 
a CISO phase, i.e. a Luther-Emery Uquid which is the univer- 
sality class of usual doped two-leg ladders. The critical fields 
at which the system changes from 4 ^ 2 or 2 ^ 4 Fermi 
points can by computed analytically (see Appendix lAli. Note 
that from Refs. I4l l33ll34i in the case of 4 Fermi points, other 
phases such as C2S1 or C2S2 appear when the difference be- 
tween the two Fermi velocities becomes sufficiently large to 
prevent runaway of some coupling constants in the RG flow. 
The large velocity difference implies that these phases are in 
the vicinity of the transition region between the CISO and the 
ClSl phase, where the Fermi velocity of the band that is emp- 
tying is going to zero. Moreover, this also implies that both 
the C2S 1 and C2S2 phases have a very small extent near the 
transition region. "^'^^-^'^ The above considerations apply for a 
system at a generic incommensurate filling. At commensu- 
rate filling, umklapp interactions can be relevant and lower 
the number of gapless modes. More specifically, at half- 
filling, an insulating phase COSO of the Mott type (resp. Band 
type) is expected for (f) < (po (resp. (j> > 0o)- It is also pos- 



sible to have an umklapp interaction inside the bonding band 
(if its Fermi wavevector equals tt/2) leading to either a COSl 
or C1S2 phase. The phase diagram of an isotropic ladder re- 
sulting from these considerations is given on Fig. [3] where the 
main feature is a reentrance of the CISO phase at high flux. 
The high flux CISO phase has a band structure very similar 
to the one of the Hubbard chain with a next-nearest hopping 
term.-'^'-'* t' for t' > t/2 and the same competing orders as 
the low flux CISO phase as we will see. This phase diagram is 
generic for t_L/i|| < 2 except that 0o < 0cWhen<^/<|| > \/2. 
For t±/t^^ > 2, the CISO phase at low flux around half-filling 
disappears. 

B. Strong-coupling limit: numerical results on the t-J model 

We now let the interactions go to the strong-coupling 
regime U ^ t where the Hubbard model (O reduces to the 
t-J model with J ~ At'^/U. In this limit, we only use isotropic 
couplings i = t|| = t± and J = J|| = J± so that the t-J 
Hamiltonian simply reads 

Ht-j - VHtV + J ^ [S, • Sj - ^n,n,] , (4) 

where Si is the spin operator and rii — c\ ^Ci_a is the elec- 
tronic density operator (leg index is omitted in and rii). V 
is the Gutzwiller projector which prevents double occupancy 
on a site. Observables are computed with DMRG for the range 
of doping 0<5 — 1 — n< 0.5. The phase diagram will be 
discussed for the special case J/t — 0.5 for which the system 
has dominant superconducting fluctuations 

1. Orbital susceptibility 

The results for the non-interacting system of Appendix lAl 
show that the orbital susceptibility plotted in Fig.fTSlchanges 
sign at the transitions from 4^2 Fermi points with sharp 
discontinuities (for Q < 5 < 0.5). It is important to note 
that the noninteracting orbital susceptibility contains contri- 
butions from all the occupied states and not just those at the 
Fermi level which control the low-energy properties. There- 
fore, such connection of the change of sign of the orbital sus- 
ceptibility and a change in the number of Fermi point is not 
obvious. Nevertheless, we propose to extend this way of prob- 
ing the phase diagram to the interacting situation. Indeed, we 
can compute the screening current j[| and its associated sus- 
ceptibility x°* as a function of the flux (p using the definitions 

and X = 

in which Eo{<j)) the ground-state energy and L is the length 
of the ladder With this definition, x°'^^{4') > corresponds 
to orbital diamagnetism. The first relation is a consequence 
of the Feynman-Hellman theorem and the second one results 
from the definition of the susceptibility as dj\i/d(f). This 
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FIG. 4: (Color online) Orbital susceptibility, charge and spin gaps 
for J/t = 0.5 and on the 5 — 0.063 line of Fig.[5] The zeros of the 
susceptibility precisely probe the different two phase transitions. 
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FIG. 5: Phase diagram of the t-J model for J/t — 0.5, for which 
the system has dominant superconducting fluctuations at zero flux, 
determined from the zeros of the susceptibility on a system with L = 
32. Results are very similar to Fig.[3]with a transition to a CIS 1 phase 
at intermediate flux and a reentrance to a C ISO phase at large flux. 



screening current can be related to the mean value of the 
current operators ji.2 along the two chains by noting that 
— O2 — ji)/2- Numerically, these quantities are com- 
puted directly from centered energy differences (to minimize 
discretization effects) 

^{<f>) = -[Eo{ct' + d(f>)-Eo{cf>-d(f>)]/{2Ld(t>) (6) 
X^'^W = [j||(0 + d<^)-j||(0-ci</.)]/(2d0), (7) 

using the conditions j|| (0) = j|| (tt) = (see AppendixlBll and 
the right and left derivatives for x"''(0) ™d x"''(7r). These 
quantity are easy to compute numerically and are found to 
have small finite size effects for J/t = 0.5. The effect of 
interactions is to smooth the discontinuities at the transitions 
but we still expect that the sign-changes of the susceptibility 
do correspond to transitions between CI SO and ClSl phases 
even in the strong-coupling limit (see Fig.|4]l. The phase dia- 
gram obtained from this ansatz is consistent with the behavior 
of other observables such as spin-spin correlation functions as 
will be seen in the next paragraphs. Thus, we can sketch on 
Fig. |5]a phase diagram similar to the one of Fig. [3]for the t-J 
model with J/t = 0.5. Compared with the weak-coupling 
phase diagram, the ClSl phase is slightly wider at low dop- 
ing but thinner for S ~ 0.5. Even in the presence of strong 
interactions, the overall shape of the phase diagram is not af- 
fected (although the precise location of the phase boundaries 
in the S, (j) plane does depend on C//t or J/t). Hence, for typ- 
ical densities 0.5 < n < 1 in the isotropic t-J model, the 
leading effect which governs the phase diagram is the change 
of the band structure under the applied flux. Lastly, one must 
note that the C 1 SO phase persists at small flux at quarter-filling 
6 = n = 0.5 in contrast to Fig. [3] This can be qualita- 
tively explained by noting that renormalization group studies 
on coupled chains have shown that interactions reduce the in- 
terchain hopping integral t± with respect to its non-interacting 
value 



2. Elementary excitations : pairing energy and spin gap 

The elementary excitations of a doped two-leg ladder are ei- 
ther the creation of a magnon, which cost is the spin gap A5, 
or the creation of two quasi-particles by breaking a Cooper 
pair, which cost defines the pairing energy Ap. Following 
Refs. i2Cl39 . we compute them numerically from the defini- 
tions 

As = SoK, S^ = l)- Eoiuh, = 0) (8) 
Ap = 2SoK-l,5^ = l/2) 

-^oK,5" = 0)-^oK-2,5^ =0) (9) 

with E^inin S^) the ground state energy of a system with rih 
holes and spin along the z axis. Since one can always 
have a 5^ = 1 state by breaking a Cooper pair, the condition 
A5 < Ap must be satisfied in the thermodynamic limit. For 
(f) — Q and J /t — 0.5, it is known-*^ that the pak-breaking 
excitation is larger than the magnon excitation (lowest triplet 
excitation). These elementary gaps as a function of the flux 
are displayed on Fig. ID We observe a decrease of the pairing 
energy with toward at the critical flux corresponding to 
the onset of the ClSl phase. In the ClSl phase, we have a 
metallic Luttinger liquid phase with zero pairing energy. The 
cancellation of the pairing energy is thus the result of a band 
emptying mechanism and should not be confused with a mag- 
netic superconducting critical field i/c2 which would corre- 
spond to a high density of vortices in a 2D (albeit anisotropic) 
superconducting material. Indeed, we do not observe a Rd 
superconducting critical field or commensurate vortex phases 
as in the model of bosonic of Ref.'24'. A situation correspond- 
ing to a true Elc2 critical field might rather be a small flux 
through an array of coupled ladder In this respect, the ap- 
proach of Ref. copmutes coiTectly the i7c2 field up to a few 
approximations. 

The spin gap increases at low magnetic field until it crosses 
Ap (see Fig. |4|i. From local hole densities (data not shown). 
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FIG. 6: (Color online) Spin correlations S{x) for fluxes in the ClSl 
phase (a) and in the CI SO phases at low and high fluxes (b) in the 
phase diagram of Fig.|5]on the S — 0.25 line. Insert: the behavior 
of the spin gap computed on the same system using Eq. {Sjl. These 
independent observables confirm the reentrance of the CI SO phase. 



the domain of hole pairs slightly shrinks which can be inter- 
preted as a reinforcement of the spin-liquid background at 
low magnetic field, in agreement with exact diagonalization 
results previously discussed.^' For larger flux but still in the 
CI SO phase, the spin gap A5 becomes identical to the pair- 
ing energy Ap and both decrease towards zero as the flux 
is increased. The energy difference Ap — A5 can be inter- 
preted as the energy of a bound state of a magnon and a hole 
pair. This magnetic resonant mode was discussed previously 
at zero flux by varying interactions-'''*^'"^' and its origin was 
related with the opening of doping controlled magnetization 
plateaus. '^'^^ Thus, the effect of adding Zeeman coupling at 
low flux {(f) < 7r/3 for J/t — 0.5) would give very similar 
results to those of Refs. Il9ll20l since the bound state survives 
to rather high magnetic flux. Finally, a small spin gap is re- 
covered at high flux (near = tt) in agreement with the weak- 
coupling limit predictions. 

To gain further insights on these excitations we have com- 
puted the spin and pair correlation functions in the ground 
state. The spin coiTelations S{x) = (S(x) • S(0)) are short- 
range in a spin-gapped phase with a coiTelation length f ~ 
I/A5 which gives complementary estimation of the evolu- 
tion of the spin gap, particularly important when the pairing 
energy is smaller than the spin gap. From Fig. |6] we find a 
similar increase of the spin gap (smaller Q at small flux and 
algebraic correlations in the ClSl phase. The spin gap at high 
flux is again recover with short-range spin coiTelations. 

The singlet operators on a rung defined by A(a;) — 
Cx,i,'\Cx,2,i — Cx,i,iCx,2,-\ gives the pairing correlations 
P{x) = (A(x)A''^(0)). While P{x) remains algebraic in the 
C 1 SO and ClSl phases, its overall magnitude follows the pair- 
ing energy and is strongly reduced in the ClSl phase. These 
correlations can be discussed from the bosonization approach 
using the conventions and definitions of AppendixIC] We find 
that the low-energy dominant term at wave-vector g = in 
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FIG. 7: (Color online) Superconducting correlations P{x) for three 
values of the flux corresponding to each phases encountered on the 
S — 0.063 line of the phase diagram of Fig. |5] The decrease ex- 
ponent is much higher in the ClSl and high-flux CI SO phases and 
the overall magnitude is strongly reduced suggesting a metallic state. 
Note that correlations display oscillations associated with a-Akp con- 
tribution. 



the low-field CI SO phase reads 

A{x) OC [{ad)^tpd,R,ai^d.L.-a + {bdfi>d.L.a1pd.,R.~a- 

a 

For instance, the intraband terms read 

i>v.R,oi^v-L-o ~ e'[ec++pflc-~^(0.++P0,-)] ^ (10) 

with j3 ± for dju. From previous results)^ we know that 
in a CI SO phase all the fields except (/)c+ are gapped, with 
(6'c-) = and (0s+) = 7r/2, (0s-) = i^jl. These terms 
are thus algebraic with a decay exponent l/(2i4rc+) which is 
the continuation of the zero-flux physics. We observe from 
Fig.|7]that /fc+ increases with the magnetic field but a precise 
evaluation is numerically difficult. 

In the ClSl phase, superconducting correlations are ex- 
pected and found to be algebraic with an exponent K'^ ^ + 
but with a much smaller amplitude, in agreement with a metal- 
lic phase. 

The physical properties of the high-flux CI SO phase are 
very similar to that of the low-field CI SO phase. Following 
the analysis by Fabrizic^^ (see Appendix ICl for notation), the 
pairing order parameter in this phase reads : 

A(a;) OC [(ai)^'(/'l,i?,(TV'l,L,-cr + {i>\f'i'\.L..cr^\,R-a 

a 

+ i^lTi'l,R,a^1,L-a + [a-lf ^l^.L ..cr^1,R-a\ ■ 

which will give fluctuations with an exponent \l(2Kc^). 
Computing the density and transverse current order parame- 
ters shows that they have a 2_ft'c+ decay exponent associated 
with the wave-vector 2{kp^\ — kF,2) = 27rn. The compet- 
ing orders are thus the same as in the low field phase studied 



in Sec. |III1 The fact that the SC signal is small and with a 
large exponent in Fig. |7] suggests that the CDW fluctuations 



o ■S « 

O ~ " CD 

dominates in this strong-coupling regime making the high flux X ■S c* I 

CI SO phase a spin-gapped metallic phase with strong trans- ^ c'-^ | 

verse current fluctuations. !c « g r-4 

* g d 

m. LOW-FIELD PROPERTIES OF THE LUTHER-EMERY ' " 



PHASE •S S 

o ^ >t 

This section discusses the properties of the CI SO Luther- ■« g .g* 

Emery phase at very low fluxes relevant to the experimental .2 g S 

accessible magnetic fields. g & 



Jn(X) — line '-■x+l,l.a'^x,l,cr <i '^x.l,a'^x+l,l,r' 



of local parallel and transverse current operators are computed ^ « S 

using the definitions =c 4^ 

1+1,1. a^x,l,cr — c ■ L. I Vx+l,l,a\ 3 9 

J\\\^) — "-IIP t-a;+i_2,CT'-2;,2,(T e L^^2,(T''a:+l,2,CTj -2 C 

•S 5^ '3 

u 



(J 

« X) 



They are related to the current operator jp along the chain p & 
via jp — j- J2x l-^) ^^^^ '■'^^ screening current (</>) . S -2 x 

We have checked that Kirchhoff 's conservation law for charge ^ 
currents is satisfied at each vertex of the lattice. At zero flux, - ^ 

no currents are present in the ladder When the magnetic field 
is applied, time-reversal symmetry is broken and local cur- 
rents have a non-zero expectation value depicted in Fig. |8]for ^ ^ 
a 2 X 32 ladder with 4 holes. First, the two hole pairs manifest 
themselves by two domains (areas with open circles). The lo- ^ ^ 
cal screening currents develop inside these domains and not at 
the edges of the ladder Clearly, in the strong-coupling limit 
where double occupancy is prohibited, the only domains in 
which electrons can take advantage of the flux are the room H c 
left by holes. This results in a periodic pattern for hole cur- o g 
rents whose length scale is exactly S^^ (— 16 in Fig. |8]l. 
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Such a length scale is different from the usual magnetic length — a ^ 

^ft/ei? governing orbits of Landau levels. A similar length § g c 
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A. Current densities and correlations S c "S ° 

'o s " ^ 

With open boundary conditions used in DMRG, we have ^ Z o ^ 

access to the local density of holes and currents. The local " i -a 

hole density reads h{x) = 1 — {n{x)) while the mean values ^ eS " * 
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scale has been found in the study of OAF phases while the 

current pattern is different from the one observed under mag- ^ S ■S 

netic field. Note that these diamagnetic currents are not re- u "3 ^ "3 ^ 

lated to the Meissner effect expected in a superconductor or in S ^ S 

a bosonic ladder.— Actually, hole pairs are delocalized on the ^ ^ Z 2 ^! 

two chains so that this pattern does not correspond to currents j o c ^ -2 

of pairs but rather to currents inside pairs. ^ c o J 2 

To study the nature of the current fluctuations, we have -2 .S ^ 

computed the transverse current correlations »j Is 

5 "o ■S c o 

j(x) = {jA^)j+m - o±(-^))(.?±(o)) , (12) 

where we have subtracted the finite local expectations (as 

one would do with density correlations). The main result of E -a c3 ^ 



o j3 M -g _g 

^ o .> g ^ 

. . 60 ^ 
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Ref. 1 1 was the absence of algebraic transverse current corre- 
lations in a CISO phase because of the strong spin fluctuations 
associated with the spin gap. Here, although the spin-gapped 
CISO phase survives at low flux, the situation is quite differ- 



ent because the chain exchange symmetry is explicitly broken 
by the magnetic field. Indeed, Eq. (fTTT i can be rewritten in the 
basis of the d, u bands, which leads in the continuum limit to 
terms with different wave-vectors: 



= it 



where the coefficients Ud/u, bd/u are defined in Appendix ICl 
We now turn to the bosonization representation of the opera- 
tors appearing in Eq. (fT3T l. Terms with the lowest wave-vector 
kd — ku contain operators of the form 



(14) 



with (T = (t,i) = ±- In the CISO phase, their correlation 
functions decay exponentially as they involve the dual fields 
6s- and (j)c- which are disordered. For the terms with the 
wave-vector kd + we find similarly that 



„i[0c++ec-+'T(0,++e,_)] 



(15) 



which also have short ranged correlation functions because 
of the presence of the disordered field 6's_ in the bosonized 
representation Eq. ( fTsl l. This is exactly the same result as 
in Ref 1 1 albeit extended to nonzero flux. Without magnetic 
field, DMRG calculationsi^ showed that the dominant wave- 
vector in the exponential signal was k'^+k'^ rather than k'^—k'^ 
probably because the bosonized expression of the correspond- 
ing Fourier component contains two strongly fluctuating field 
for the latter term (see Eq. ( fT4l ) , but only one (see Eq. ( fT5] )) for 
the former term. In the presence of a magnetic field, we see 
that the terms in Eq. (T3[ corresponding to the wavevectors 
2ku and 2kd have a magnitude 



(bp 



at small flux 0, i.e. they exactly cancel for = but are 
present once the magnetic field is turned on. These new term 
are allowed by the symmetry reduction induced by the mag- 
netic field. Furthermore, they have the bosonized expression: 



(16) 



in which the spin contribution has long range order, but the 
charge contribution contains the dual field 0c- leading to ex- 
ponential decay of the associated correlation function. There- 
fore, all the "2fci?" contributions in the transverse current cor- 
relation functions display an exponential decay. 

In the expression ( fT3] l, the higher "Akp" harmonics of the 
current are not taken into account. The reason for this is that 
in the full Hilbert space, the Akp component of jj^ is simply 



proportional to c\f^^C-kF- However, in bosonization, a mo- 
mentum cutoff is introduced, and the high energy states which 
involve the creation or annihilation of fermions with momen- 
tum farther from the Fermi momenta than the cutoff are elim- 
inated. Thus, the expression of the current at lowest order in 
the interaction U/tin Eq. ( fT3] l cannot contain any Akp contri- 
butions. However, virtual processes in which a fermion is cre- 
ated and annihilated far from the Fermi points give contribu- 
tions of higher order in U/t to the Akp components that only 
involve fermion operators close to the Fermi points. Such con- 
tributions can be derived in perturbation theory following the 
approach in Refs i42i43l In the case of the transverse current. 



an interaction of the form Uci 



C3k, 

.t 

■k, 



yields a pertur- 



bative contribution proportional to U/t c\.^cl.^C-kpC-kp to 
the Ak F component which involve only operators belonging to 
the low energy subspace and thus cannot be neglected. Such 
corrections can be viewed as a pair hopping or a correlated 
hopping between the chains. Thus, we expect to find ?i Akp 
contribution to the transverse current of the form: 



(17) 




J/t = Q.5 and 8= 1/16 



5x10 



-5x10 



10 20 30 40 50 60 
X 



FIG. 10: (a) Transverse current correlations become algebraic in the 
CISO phase once the magnetic field is turned on (same parameters as 
in Fig. [8}. (b) Demodulation of the signal enables to extract clearly 
the wave-vector 2n5 of the correlations. 



TABLE I: Summary of the bosonization result for operators in the 
low field CI SO phase (see Appendix [Cl for notations). We have 
2kF = Trn. If not short-range ("exp." notation), we give the decay 
exponent of the associated correlations. Numerically, {n{x)n{0)) 
and {j± {x)j± (0)) are algebraic because they pick up the Akp terms. 



Operator 


in the CISO phase 


exponent 


wave-vector 




exp. 


2kF 


A{x) 


l/(2^c+) 





n2kp(x) 


exp. 


2kF 


n4kp{x) 


2Kc+ 


4kF 


j±,2kp (x) 


exp. 


2kF 




2K,+ 


4kF 



associated with the wave-vector 2{k'^ + fc") = 27r(l — S). 
The {j^^kp {x)j±,4:kF (0)) correlations have a power-law de- 
cay with exponent 2Kc+. This result is very similar to the 
CDW fluctuations associated with the correlations {n{x)n{Q)) 
which, while being short-range at 2kp, also possess a Akp 
power-law decay. ^ Note that these CDW correlations con- 
tain terms analogous to (fT4l [TSl \T6[ but with different pref- 
actors. This is in good agreement with numerical results for 
which we found algebraic correlations with wave-vector 2tt6. 
The wave-length of the coiTelations is again associated to the 
length scale S^^ of the local hole and transverse current pat- 
terns. We found numerically a larger Luttinger exponent from 
the current correlations Kc+ ~ 1 while superconducting cor- 
relations rather give Kc+ ^0.6 for the same parameters. This 
difference, also found for charge correlations,^ could be at- 
tributed to the need for virtual high energy processes to create 
the Akp coiTelations, leading to somehow low and noisy sig- 
nals. The behavior of the bosonized operators is summarized 
in table U 
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0.25 0.5 



5 

FIG. 1 1 : (Color online) Zero field susceptibility xo as a function of 
doping S and interactions parameters J/t for a ladder with L = 48. 
A strong reduction at low J/t and commensurabilities 5 = 1/4 and 
1/2 are clearly visible. The computed susceptibility is nearly zero at 
S — 1/4 but subject ot finite size effects (see Appendix iBt. On the 
contrary, the 5 = 1/2 phase has a finite susceptibility. Insert: the 
derivative dxo/d6\s=o as a function of J/t. 



creases continuously as J/t is lowered (data not shown, for 
larger J/ 1, the convergence is better) so that we are confident 
that the observation is not an artifact of the finite size effects. 
For 6 = 1/2, a discontinuity of the slope is found but the 
susceptibility remains finite in this phase (there, the finite size 
effects are smaller). The occurrence of insulating CDWs was 
previously studied^ in this part of the (5, J/t) phase diagram 
of the t-J model. However, only the (5 = 1/4 CDW phase was 
discussed and Fig. [TT] suggests that the 6=1/2 phase is of a 
different nature. 



B. Zero-field susceptibility and tlie commensurate phases 

Since only very small flux per plaquette cf) can be achieved 
experimentally, we now focus on the limit of the or- 

bital susceptibility xq = X°'''(0) which is calculated numeri- 
cally from Eq. (|5]l. For the non-interacting system, this quan- 
tity is finite and positive at half-filling and then increases with 
doping (see Fig. [16] of Appendix [All. Once interactions are 
turned on, this susceptibility is zero at half-filling because of 
the Mott insulating state (see Fig. [TT] and Appendix [A] for a 
general discussion of the susceptibility). When doped, the 
system acquires a susceptibility roughly proportional to den- 
sity of charge carriers with xo ^ 5- The proportionality coef- 
ficient decreases with J/t which is reminiscent of the fact that 
large Js reduce the mobility of holes (see insert of Fig. [TTT l. 
Compared with the non-interacting result, the susceptibility is 
thus strongly reduced by the interaction. When the ratio J/ 1 
is lowered, a strong reduction of xo is clearly visible for the 
hole commensurate doping 6 = 1 /4 up to finite size effects 
discussed in Appendix IB] This drop of the susceptibility in- 
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FIG. 12: (Color online) Two-particle charge gap as a function of 
doping showing the emergence of insulating phases at commensura- 
bilities 5=1/2 and 5 = 1/4 at low J/t. The smallest flux we have 
access to is enough to destroy the 5 = 1/4 CDW phase (see text 
for details), while the S = 0.5 BOW phase is more stable. Data are 
computed with an extra J = +0.3 at the two extremal rungs. 
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To study more precisely the occurrence of these phases, we 
have computed the two-particle charge gap 

A2p = En{nh + 2) + Enim, - 2) - 2Eo{ni,) (18) 

as a function of doping and interaction at zero magnetic field 
and for the lowest flux Att/L we have access to. A system 
with pairs always has a finite one-particle charge gap (or pair- 
ing energy), but is insulating if the two-particle charge gap is 
also finite. The results on a ladder of finite length L = 48 
are given on Fig. [12] On the one hand, a strong disconti- 
nuity at (5 = 1/2 is found even for rather large J/t and the 
charge gap of this phase survives to nonzero flux (away from 
the commensurability, is much smaller but finite because 
of the finite length of the system). On the other hand, the 
(5=1/4 discontinuity is only found at small J/ 1 and is de- 
stroyed for the lowest flux we can use. No other discontinuity 
of the two-particle charge gap is found for the range of doping 
< S < 0.6. The system with S — 1/4 has edge effects and 
we have added an extra J± = +0.3 on the two extremal rungs 
to control the spinons at the edges as it was done in Ref 6. This 
phase is thus difficult to study under magnetic flux but it was 
studied previously and it was proposed to be a four-fold de- 
generate CDW phase with pairing and a small spin gap on the 
basis of the behavior of the Friedel oscillations. These gaps 
were found to be numerically very small. ^ The observed sen- 
sitivity to the flux is consistent with small gaps. Indeed, such 
a four-fold degenerate phase is difficult to stabilize. Qualita- 
tively, if pairs of holes are well-formed on rungs, it is hard to 
generate an effective long range repulsion between these pairs 
to stabilize a crystal of hole pairs. On the contrary, if pairs are 
spread over a few rungs, they can repel each other more easily 
but will have a smaller spin gap and pairing energy. This latest 
picture of a pair of holes delocalized on a plaquette every two 
plaquettes seems to be more suited to describe this phase. 

In the insulating phase with S = 1/2 (quarter-filling), 
instead of the pronounced Friedel oscillations obtained for 
6 = 1/4, a uniform electronic density (see Fig.|9]l is found. 
However, if one computes the bond order parameters t^, (x) 
along the bonds at zero magnetic field by using the definitions 

^1,11(3;) = ^||(4+l,l,<TCa;a,^ +4,l^^C2,+l,l,^) 
^2,||(a;) = i||(4+l,2,<TCx,2,o- +4_2,<tC2.+1,2,<t) 
tjLix) = i_L(4^i^^C:r^2,a +4,2,^C:r,l.<T) , 

one finds strong oscillations with a period of two lattice sites 
for the II bonds making this phase an insulating bond or- 
der wave (BOW) phase (see Fig. [TsT l. This is confirmed 
by the current pattern under magnetic field found in Fig. |9] 
which has well-defined orbits around plaquettes but small 
currents between plaquettes. The local transverse current is 
staggered while the transerve bond density wave order pa- 
rameter is uniform. Such local orbits allow a finite orbital 
susceptibility even if the system remains insulating because 
0±(a;)j±(0)) ^ (j±(l)j±(0)) in Eq. (lD2j which is simply 
the local response on a plaquette. 




6 12 18 24 

X 

FIG. 13: (Color online) Local kinetic bonds t^{x) in the 5=1/2 
bond order wave (BOW) phase computed at zero magnetic field. The 
parallel bond orders are strongly oscillating at wave- vector n while 
the transverse bond and the transverse kinetic bond and electronic 
densities are uniform (see Fig.|9j. 



IV. CONCLUSION 

We have studied the effect of a magnetic flux through a 
doped two-leg ladder by means of bosonization and DMRG 
calculations. As a function of the flux, a rich phase diagram is 
observed with an intermediate Luttinger liquid phase and the 
reentrance of the Luther-Emery phase at high flux. Both in the 
weak- and in the strong-coupling limits, the phase diagram is 
governed by the evolution of the band structure. Focusing on 
the small field physics of the Luther-Emery phase, we observe 
that local currents develop in the ladder inside the hole pairs 
regions. Their typical length-scale S^^ is controlled by hole 
doping S. The transverse current correlations also develop 
as soon as the magnetic field is turned on with an algebraic 
behavior contrary to what was found without magnetic field. 
Lastly, we have computed numerically the zero-field suscepti- 
bility of the system as a function of the interaction parameter 
J/t and of hole doping. We found that insulating commen- 
surate phases at low J/t exists only at dopings 6=1/4 and 
1/2 and that the two phases have different responses under 
magnetic field. The contribution of the conduction electrons 
to the orbital susceptibility might thus be useful to probe these 
phases experimentally. Results on the 6 = 1/4 phase are con- 
sistent with a four-fold degenerate ground state with a small 
pairing and spin-gap, making it very sensitive to the flux. On 
the contrary, the 6=1/2 phase appears to be a robust bond or- 
der wave phase with a two-fold degenerate ground state. De- 
spite its insulating nature, this phase has a finite susceptibility 
due to local orbits of electrons around plaquettes. 

The ladder compound Sri4_2:Caa;Cu2404i (SCCO) was 
the first non-square cuprate compounds showing super- 
conductivity under high pressure."*"* The presence of a 
spin gap in its superconducting phase has been addressed 
experimentall y'*^i"*^i"*^'" *^ but no consensus has risen on the 
actual nature of superconductivity in this material. An- 
other exciting feature of SCCO is the occurrence of charge 
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density waves at ambient pressure i'*^i^"'^^'^^-^-^'^^-^^-^^'^^i^^ 
Experiments^^ ?^ have suggested that CDW could appear at 
wave-vectors 1/3 and 1/5 which was discussed theoretically 
using a multiband charge transfer model solved by Hartree- 
Fock approximation.^'' Here, we have showed that the t-J 
model on a single ladder only displays 1/4 and 1/2 commensu- 
rabilities, as proposed in Ref. 6, and that orbital susceptibility 
could help to understand the nature of these commensurate 
phases. 

We would like to point out that an interesting realization of 
quasi-one dimensional systems in which magnetic flux can af- 
fect the band structure is provided by carbon nanotubes.^" Fol- 
lowing the theoretical prediction,*" experiments on multiwall 
nanotubes (where notable fluxes can be achieved due to the 
large diameter of the outmost shell) have shown that the band 
structure of these systems was indeed sensitive to magnetic 
fluxesi^ii^2iS In the case of gapped zig-zag single wall nan- 
otubes, although the experimentally accessible fluxes through 
the tubes were small, an effect on the conductance oscillations 
in the Fabry-Perot regime could nevertheless be evidenced as 
a result of the lifting of the degeneracy between two subbands. 
As there exists some evidence for strong electronic correla- 
tions in carbon nanotubes^i^Si^^i^LS and as carbon nanotubes 
possess some analogies with ladders^^^ an interesting exten- 
sion of the theoretical results developed in the present paper 
would be to study quasi-one dimensional models mimicking 
more closely carbon nanotubes. It would be particularly in- 
teresting to compute the behavior of the Luttinger exponent 
controlling the Zero Bias Anomaly as a function of the ap- 
plied field. 
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Hamiltonian by taking the Fourier transform, which gives the 
energies 

£d/u{k,(l)) = -2t|| |cosfccos^ ± y^sin^ fcsin^ ^ (f 



(Al) 

The basis transformation can be written using coherence fac- 
tors afe,6fc > 




with 



at bk \ I Ck,d 



sin k sin ^ 
sin^fcsin^ f + (f] 

sin k sin $ 



sin-'fcsin'^ f + (f) 



(A2) 



(A3) 



(A4) 



We have a-k = bk and the factors depend on the wave-vector 
k while at zero flux ak = bk = l/\/2. For any finite flux and 
fc > we have < bk with 



1,2 2 

h - ak 



sin k sin ^ 



(A5) 



Below are a few considerations on the band structure ( lAll ) 
which are summarized in Fig. ID 

(i) The condition to have a band gap is that maxEd — 
5(i(7r, 0) < mrnSu = ^uiO, (f>) which gives j = cos |. This 
condition can be reformulated as 4> > (po with 



sin — = W 1 



(f 



(A6) 



(ii) The condition to have a double well in the £d band 
comes from the sign of the second derivative at fc = and 



APPENDIX A: FORMULAS FOR THE NON-INTERACTING 
SYSTEM 



Part of these results were first given in Refs.mH We re- 
produce them for clarity and notation conventions (which are 
different) and extend them when necessary. In what follows 

1. Band structure 

At zero flux, the interchain coupling lifts the degeneracy be- 
tween chains giving birth to a bonding band ky = and anti- 
bonding band ky = tt. The flux breaks the reflection symme- 
try between chains and couples these and tt modes. We call 
down and up (with labels d/u) the two bands in presence of 
a flux. It is straightforward to diagonalize the non-interacting 




FIG. 14: Depending on the filling, the number of Fermi points can be 
either 2 or 4 (sketched on the bands of the B phase of Fig.[2]l. Note 
that at low and high filling, the two Fermi velocities have opposite 
signs while at intermediate filling, they have the same sign. 



11 



is ^ cos ^ = sin ^. This condition can be reformulated as 



> 6c with 



sm ■ 



1/2 



(A7) 



In this case, the wave-vectors corresponding to the minimum 
of the down band and the maximum of the up band read 



k^ini'f') = arcsin 



(f 



cot^ 



(A8) 
(A9) 



Working at fixed electronic density n, we have to relate the 
Fermi points directly to n using Eq. dAll ) and the Luttinger 
sum rule (see Fig. [14] for sketch of all possible situations). If 
the system has only two Fermi points, then we have either 
kp = Trn or hp — Ti{n — 1), where both relations do not 
depend on the flux. When there are four Fermi points and 
limiting the discussion to n < 1, we have either that: 
(i) the bands are overlapping: 



kl 



k% 



nn 



sm 



h,p 



■ 2 d> 

sm f 



• 2 V 

sm — 
2 



1/2 



which appears with a finite value. 

(iii) The two curves intersection is a = \/2 for which we 
have 00/c = n/2. \f a < ac we have 0c < <t>Q, else (pa < (f>c. 

(iv) The condition to empty the u band is pt ~ a — 2 cos ^ 
and is the same as the condition for which the d band is com- 
pletely filled. 

(v) We can show that only situations with two or four Fermi 
points can occur: for monotonous bands (cj) < (j>c), this is 
obvious. For non-monotonous dispersion ((/> > (f>c), if (f) > (po, 
the two bands are not overlapping so we have either two or 
four Fermi points. For any flux (p G [<Pc, '/'o]^ we can convince 
ourselves that since £"{k, (f) > £'^{k, (f) (for any k) and since 
the up band has a unique maximum (at fcmax(</*)) while the 
down band has a unique minimum (at k'^^J^cp)), it necessarily 
implies that only two or four Fermi points are allowed (see 
Fig. [Hi. 



2. Filling the bands : finding Fermi points 



Fermi points k'^l^ are deduced from their relation to the 
chemical potential /i from Eq. ( lAll l. If needed, this equation 
can be inverted into 



cos kp^"{fj,, (p) 



- t: — cos — 
2U 2 



± 



Sd/v 



2^11 



sm 



(p fOi\'^ 

2 + 12) ■ 
(AlO) 



(ii) there is a non-monotonous dispersion, for the down band: 



1 ' F.2 



(P^cos^l 
sin2(^) - sin^ 



1/2 



These equations reproduce the correct result in the (p ~ and 
a = limits. It is also straightforward to compute the 4 crit- 
ical densities at which the number of Fermi points changes 
from 2 to 4 and 4 to 2 (see Fig.|3]l. By arranging them accord- 
ing to fici < < < rici, we have 

rici — arccos(cos0 + acos(0/2))/7r if e [i/ic , tt] 

nc2 — arccos(cos0 — acos(0/2))/7r if (/> g [0, (/)o] 

= 2 - if e [0, (pa] 

nc4 = 2 - rici if (/> e [0c, tt] (A12) 

These equations can be inverted to give the critical flux 0^*^^ 
at which the transitions from 2 <-> 4 Fermi points occur: 



A2^4 



±- 



4 cos- 



/'Trn\ 



(A13) 



this last expression is useful to check that the change of sign 
of the susceptibility is associated with these transitions. 



3. Orbital susceptibility 



This latest is useful when working at fixed fi. Note that, de- 
pending on (p and /i, the above equation can have two roots 
labelled by s^/u = ±1 for each sector d,u (see Fig. [T4b. If 
there are two Fermi points in the same band p, we use the 
notation fc^ The Fermi velocities can be evaluated from 



v'p^ik, (p) = 2i|| sin A; ^ cos =F 



cos k sin ^ 



sin-^fcsin^l + d) 



(All) 

Note that one of the two Fermi velocities is negative in case 
of non-monotonous bands. 



Knowing the location of the Fermi points , we can com- 
pute quantities integrated over the bands such as the total 
energy Eo{(p), the screening current j\\{(f>) and the associ- 
ated orbital susceptibility from Eqs. (|5]l. One contribution to 
this current from electrons of momentum k is ^^^{k, <p) — 
9£d/u{k, (p)/d(p, which reads: 



■ 11 ^ cosfcsm- =F 2' 



sin k sin ( 



sin^ k sin^ ■ 



(A14) 
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FIG. 15: Susceptibility as a probe to the number of Fermi points for 
0.5 < 1 — 5 < 1.0 and ~ - 0^"'* indicate the transitions from 
4 to 2 and 2 to 4 Fermi points. 




-0.5 - 



0.5 1 

n 

FIG. 16: Zero field orbital susceptibility of the non-interacting sys- 
tem for t± = t|| from Eqs. ( IA16IA171 ). The singularity at n = 0.5 
marks the transition from 2 to 4 Fermi points. When n < 0.5, the 
susceptibility changes sign for n ~ 0.3975 . . . while the system al- 
ways has two Fermi points. 



and similarly, the corresponding contribution to the suscepti- 
bility reads 



Xd/n(fc- (/)) = - ^ COS fc COS I T 



sin^ k cos < 



1/2 



±- 



• 4 T -2 

sm K sm I 



sin^fcsin^f + (i; 



3/2 



(A15) 



It is important to remark that additional contributions come 
from the derivatives dkP{(j)) /d4> and d'^k'P{(t)) / d(tP' because, 
in the case of four Fermi points, the location of the Fermi 
points depends on the flux. Thus, Eqs. (|5]i can be computed 
either analytically or numerically. A typical plot of the inte- 
grated x°*{4>) is given on Fig. [TS] which shows that its dis- 
continuities are associated with the 2^4 Fermi points tran- 
sitions of Eqs. ( |A13l l. We have checked that the latest result 
is valid only for 0.5 < n = 1 — 6 < 1.0 when a — 1 (see 
zero-field susceptibility below). 

Lastly, the behavior of the zero-field susceptibility Xo^ can 
be computed as a function of the electronic density n. With a 
factor two for the spins, we have: 
if a > 1 — cos nn (two Fermi points). 



Xo\n) 



, 1 \ ' 

sin TTO I 1 H cos 7rn n 

a I a 



(A16) 



else, if a < 1 — cos 7m (four Fermi points), we have 



sm 



cos 7rn 



sin''(7rn/2) 



2 fa 

- — arcsm — -; — -— 

a \ 2 sin(7rn/2) 



(A17) 



The curve for a = 1 is plotted on Fig. [161 which shows that 
for n < 0.5, even when there are only two Fermi points, the 
susceptibility can be either positive or negative. 



APPENDIX B: FLUX QUANTIZATION AND FINITE SIZE 
EFFECTS 



In this section, we discuss the quantization of the flux on a 
finite size ladder First, Hamiltonian ([T]) clearly gives E{(f)) = 
E{2tt — (p) so that we can restrict ourselves to the window <j) E 
[0, tt]. Similarly, from Eqs. (|5]l we have j{(p) = — j(27r — 
which implies that j{0) = j{TT) = 0. What is quantization 
of the flux on a finite system ? Using periodic boundary 
conditions with the gauge of Fig. [T] the integrated flux along a 
leg is ±0/2 X i so that there is no remanent flux through the 
cylinder hole of the periodic ladder if 



mAn/L 



(Bl) 



with m an integer Actually, this can also be simply under- 
stood from momentum quantization k — 2TTm/L and looking 

at the dispersions —2t|| cos(fc±(/)/2) on each leg when tx = 0. 
This quantization can be checked numerically with exact diag- 
onalization. Furthermore, another possible gauge which gives 
the same flux per plaquette is to take (t)±{x) — 4>x, with 0^ 
the flux along a rung, and no flux along the legs. We can show 
that the two gauges are strictly equivalent on a finite system 
with periodic boundary conditions only if (IB II ) is satisfied. 

With DMRG, we are using open boundary conditions for 
which we expect similar effects due to momentum quantiza- 
tion. For most quantities and parameters J/ 1, the points with 
(f) = 2TTm/L interpolates nicely with the ones using dBlb so 
we can relax the constraint (see Fig.|4]for instance). However, 
when taking the derivatives such as in Q at low J/t where 
more finite size effects are present, one must strictly obey (IB 11 1 



13 



Eo((^) 
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o J/t = 


0.25 d(|) 


= 2jt/L 


D J/t = 
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0.40 d(|) 


= 2jt/L 


. J/t = 


0.40 d(|) 


= 4jt/L 
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0.25 
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FIG. 17: (Color online) Finite size effects on the calculation of the 
zero-field susceptibility at low J/t. (a) energies as a function of flux, 
(b) xo as a function of the doping S from Eq. (B21 for two different 



to have correct estimates (this is necessary in Fig. [TT| for in- 
stance). On Fig. [17] because we know that j|| (0) = 0, the zero 
field susceptibility is 



(B2) 



Taking dcj) = 2tt/L gives a negative susceptibility while d(t) — 
Att/L gives Xo ~ 0. Finite size effects are smaller away from 
the 6 = 0.25 commensurability and also at larger J/t and . 

The DMRG simulations were performed using the standard 
finite system algorithm on systems ranging from L = 32 to 
L = 64, with minor modifications to treat complex wavefunc- 
tions. The density matrix at each step was the sum of the den- 
sity matrices constructed using the real and imaginary parts 
separately. Typically, we kept m = 1400 states per block, 
giving a typical truncation error of 10^^. The presence of the 
magnetic field did not increase the truncation error notably. 
Correlation functions are computed by averaging two-point 
correlations of equal distance. Correlations with one point 
too close to one of the edges are removed. Even if there is 
no translation symmetry with open boundary conditions, this 
method gives comparable results with the one using one point 
fixed at the middle (but this method gives access to a larger 
distance x). 



APPENDIX C: BOSONIZATION CONVENTIONS 

We use the same conventions as in Ref. [t^ The bosoniza- 
tion procedure starts from the linearization of the band dis- 
persion in the vicinity the Fermi points. When there are four 
Fermi points, two of then in the up band and the two other 
in the down band down bands (corresponding to the low-field 
CI SO phase), we use: 




with implicit spin index if not explicitly required. We denote 
by ipR/L.d/u the bosonized right and left movers inside each 
bands. Note that we have different Fermi levels kp^d ^ kd 
kp^u = ku- From ( IA3b and ( IA4I ). we deduce that a-k^/^ = 
= ^d/u and h-ka,^ = at^,^ = aa/u- The bosonized 
version of the local fermion operators depends on how many 
Fermi points we have and which bands are filled. If we have 
four Fermi points and overlapping bands, we use 

C2(a;)/V^ ^ bdC^'^'^'i^dM^) + ade~''''"'ipd.L{^) 

where the right and left moving Fermi fields have the 
bosonized representation 



V': 



p,r,cr 



J<Lr<pr,p,a 



with a. a cutoff (not t_L/<||) and r = R/L, p = d/u and 
'^'R/L = T1- Vp a are Klein factors than ensure anticommu- 
tation of fermion operators having different spin or band in- 
dex. The fields (l)r,p,cr are chiral boson fields. The non-chiral 
bosons fields are defined by: 

" [4>L,p,a + 0i?,p,o-]/2, 



f^p,(T 

^p.cr 



i>L,p,cr — (l)R.p.a]/2, 



(CI) 
(C2) 



and they satisfy commutation relations [(j>p^a-{x),9p'^a'{x')] — 
i5ppi5cr.cr'5{x — x'). As usual in the framework of two coupled 
chains, we also introduce the following combinations of the (f) 
and 9 fields: the charge and spin modes in each bands p are 



0s,p = [^p,T - </'pa]/^/2 



(C3) 
(C4) 



and similar transformations for the 0. And lastly, the ± com- 
binations 



c/s,± 



(C5) 



The Luttinger parameters associated with these bosons are 
Kc± for the charge sectors and Ks± for the spin sectors. 

In the case of two Fermi points (intermediate flux ClSl 
phase) and n < 1, only the down band is filled and we can use 
the results of a single chain but using 

ci(a;)/V^ ^ ade^'^'^'^i^dAx) +bde''''''''ipdAx) 
C2{x)/^/^ bde'^''''^dM{x)+ade-'''-'''il)d,L{x). 

We simply denote by Kc and Ks die Luttinger parameters 
corresponding to the charge and spin modes. 

In the high-flux CI SO phase, four points are present in the 
down band. With the notation fcF,i,d = fci 7^ fcF,2,d = ^2, we 
have after linearizing the band structure: 

ci{x)ls/^ aie^^^^'^iAx) + 6ie-''=i^Vi,L(a;) 

+ a2e''^^V'2,L(x) + fe2e-^'=^"V2,i?(x) 
C2(x)/yS ^ 6ie^'=i-Vi,i?(a;) + a^e-'^^'^^^Ax) 

a2e-*'=^"V2,i?,(x) , 



+ 62e^'^"V'2,L(x) 
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and similar expressions for the Fermi operators ipp,r{x)- 

APPENDIX D: DIAMAGNETIC SUSCEPTIBILITY 

Let us consider the Hamiltonian (O or (01 in the limit ( 
0. By expanding to second order, we have: 



where ((;))o represents the retarded response function and ()o 
is the expectation value in the ground state without magnetic 
field. In the absence of interchain hopping, the cross response 
function (Oi;.?2))o would vanish and (j||)o would be simply 
the sum of Drude weights of each chain. The expression of 
(j|[)o can be rearranged by noting that: 



H = 7^(0 = 0) - - j2) - y i(A'i + K2), 

where ji 2 are the densities of current operators along chains 
1 and 2 and Ki,2 represent the densities of kinetic energy in 
chains 1 and 2. Note that these operators are taken at = 
0. We obtain the density of screening current operator j|| = 

-^dn/dc/j as: 

Using linear response theory, we obtain the expectation value 
of this current in this limit as: 



J\\W = 0||)0 = 5 [(((jl - J2); (jl - J2)))0 + {Kl + K2)o] 



(Dl) 



as a consequence of Kirchhoff 's law. So we have that: 



Xo - / dx{j^{x)j^{0))o + -{K, + if2>o . (D2) 



In the case of negligible transverse current correlations, this 
term reduces to the expectation value of the kinetic energy. In 
an insulator, this yields a vanishing diamagnetic susceptibility. 
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